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Abstract 

Light (M < 20 MeV) dark-matter particles freeze out after neutrino decoupling. If the dark-matter 
particle couples to a neutrino or an electromagnetic plasma, the late time entropy production from 
dark-matter annihilation can change the neutrino-to-photon temperature ratio, and equally the ef¬ 
fective number of neutrinos A^eff- We study the non-equilibrium effects of dark-matter annihilation 
on the A^eff and the effects by using a thermal equilibrium approximation. Both results are con¬ 
strained with Planck observations. We demonstrate that the lower bounds of the dark-matter mass 
and the possibilities of the existence of additional radiation particles are more strongly constrained 
for dark-matter annihilation process in non-equilibrium. 
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I. INTRODUCTION 


Photons and nentrinos are the lightest particles in the Standard Model (SM) and give the 
radiation energy density at late times in the early nniverse. The SM nentrino species con- 
tribntes three degrees of freedom becanse there are exactly three nentrino mass eigenstates 
(z/i, i> 2 , vz) in combinations of the three flavor eigenstates {ve, Vp,, Vr) of the weak interaction. 
The weak interactions that keep nentrinos in thermal contact with the electromagnetic 
plasma become ineffective aronnd a second after the Big Bang. Nentrinos deconple at a 
temperatnre on the order of 2 — 3 MeV before pairs annihilate and, thns, do not share 
in the entropy transfer from pairs. This canses the nentrino temperatnre to be less than 
the photon temperatnre later. However, nentrino deconpling was not qnite complete when 
e’''e“ annihilation began, so some of the energy and the entropy of photons conld transfer to 
nentrinos. If the dark radiation density is parameterized in terms of the nnmber of effective 
nentrino species with the canonical nentrino-to-photon temperatnre ratio, Nes increases 


to slightly more than the three nentrino species, leading to = 3.046 


flial. 


Becanse the 


nnmber of effective nentrino species is precisely predicted in the SM, this can give a ro- 
bnst constraint to any nonstandard physics. For example, new relativistic particles, snch as 
the light sterile nentrino or the Goldstone boson jd] which has a deconpling temperatnre 
less than 100 MeV, arise in many extensions of the SM, and their existence will contribnte to 
the dark radiation energy density. This scenario is, however, strongly exclnded at over the 


3ct level in the 
at a later time 


atest Planck analysis 

Ml- 


nnless photons or electrons (positrons) are heated 


According to a recent analysis of the cosmic microwave backgronnd (CMB) temper¬ 


atnre anisotropy by the Planck satellite j^, Nes was fonnd to be 3.15 ± 0.23 (lex), 
consistent with the SM prediction. We shonld recognize that the Hnbble constant 
(hfo = 67.8 ± 0.9 kms“^Mpc“^) inferred by Planck is in tension at abont 2.4cr with the 
direct measnrement of Hq = 73.8 ± 2.4 kms“^Mpc“^ by the HST lOj]: larger valnes of the 
Hnbble constant prefer larger valnes of Wg. Aeff was not strongly exclnded to abont the 
2 ct npper limit in the Planck analysis. Additionally, Aeff can be inferred from big bang 
nncleosynthesis (BBN) considerations [n| at times earlier than recombination becanse the 
theoretical expectations for the primordial abnndances of light elements depend on Wg. 
Recently, two gronps annonneed different resnlts for Wg determined from an analysis of ^He 
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13|. One group 
15| found a larger 


abundance measurements in combination with the D abundance 
obtained a result consistent with the Planck observation, but the other 
value of Ngfi ~ 3.58. This probe does not have the same resolving power as the Planck 
satellite. We will use Planck results to constrain the dark radiation in this work. 


Recently, weakly interacting massive particles (WIMPs) with sub-GeV masses 16l-l21|. 
referred to as light dark matter (DM), have received some attention because the existence of 
WIMPs with mass less than 20 MeV can modify the early universe energy’s and entropy den¬ 
sities. If DM particles couple to the SM particles (neutrinos, or photons and pairs) and 
are sufficiently light (M < 20 MeV) to annihilate after neutrino decoupling, the primordial 
plasma will be heated by DM annihilation. This will affect the neutrino-to-photon tempera¬ 
ture ratio and so might be able to explain possible differences of Wfr from the SM prediction 
or to avoid strong experimental constraints on the existence of additional radiation parti¬ 


cles. This scenario was studied in the equilibrium version for neutrino heating 0,y,i22i 
and photon heating 6-^. The equilibrium version is, however, a rough approximation be¬ 
cause DM particles are nonrelativistic at freeze-out^. If the DM particles are relativistic 
at DM decoupling, the equilibrium method must be a good approximation because they 


will decouple at equilibrium concentrations. A more accurate description of the freeze-out 
process should be considered when DM particles are nonrelativistic at freeze-out. Here, 
the Boltzmann equation is applied to the time-evolution of the DM number in a spatially 
homogeneous and isotropic universe. We treat an adiabatic expansion of the universe so 
that the total entropy stays constant and the second law of thermodynamics can be applied 
to the entropy (temperature) evolution of the produced relativistic particles. Because the 
DM mass (M) determines the dark radiation energy, it will be the parameter that we will 
constrain. We start to study the DM number evolution for an expanding universe and dark 
radiation (We) in the equilibrium approximation. Then, we study the out-of-equilibrium 
light-particle production, its entropy (temperature) evolution and its effect on the dark ra¬ 
diation. The possibility of the existence of new light species (equivalent neutrinos) is also 


^ We distinguish terminology, “decoupling” and “freeze-out”, in this paper. “Decoupling” will be used in 
the case that (DM) particles are completely non-interacting at some point, and “freeze-out” is for chemical 
decoupling. Notice this does not mean that DM dumps energy in the neutrino or the electromagnetic 
plasma instantaneously. The evolving comoving number density of DM particles has a sizable deviation 
from its equilibrium prediction around freeze-out (see Fig. 1). 
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investigated. 


II. THEORETICAL DETAILS AND NUMERICAL RESULTS 


The dark radiation energy density of the universe is parameterized in terms of the 
energy density of photons and the effective number of neutrinos iVeff with the neutrino- 
to-photon temperature ratio of the SM given by 
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Pe 


= Iiv. 


eff 


( 1 ) 


Pi ^ \-‘7/ SM 

The factor 7/8 is due to the effect of Fermi-Dirac statistics on the energy density. The exact 
formula for TVeff depends on the scenario (model). Because the temperature will change in 
our scenario, iVeff can be expressed by 


T 

-L 'll 


( 2 ) 


'-1 / \ 7 / SM 

We have assumed that DM particles are nonrelativistic at the time we consider, for example, 
the time of DM freeze-out or recombination. The BBN imposes limits on iVeff at photon 
temperatures around 1 — 0.1 MeV, and any additional dark radiation particle is unfavor¬ 
able from BBN considerations though a small possibility still exists. If the DM particle is 
relativistic in the BBN era, it becomes a dark radiation particle. We investigate the effects 
of DM annihilation on iVeff for the DM mass range of 0.1 — 20 MeV. The ratio of the neu¬ 
trino temperature to the photon temperature can be determined by entropy conservation 
because the total entropy stays constant in an adiabatic expansion of the universe. After 
neutrino decoupling, the primeval plasma will consist of two decoupled components, the 
electromagnetic component and three neutrino ones. The entropies of the neutrino and the 
electromagnetic plasmas are separately conserved, and this must serve as an efficient tool 
for the study of dark radiation here. 

Two independent thermal baths exist after neutrino decoupling, and we will consider the 
DM interaction (annihilation) in each thermal bath. Here, the DM particle always interacts 
with the plasma in thermal bath a, and thermal bath b is not relevant to the DM interaction, 
he., a, 6 = z/ or 7 , but a ^ h. 
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A. Dark-matter Number Evolution 


DM particles are in thermal contact with the rest of the cosmic plasma at high temper¬ 
atures, but they will experience the freeze-out at a critical temperature. In this case, we 
should consider the Boltzmann equation for the evolution of DM number. Because the DM 
interacts with one of the plasmas, we express the comoving number density about the tem¬ 
perature of the plasma in thermal bath b. This notation is very useful because one plasma 
can always be in the thermal equilibrium. If the DM interacts with the plasma in thermal 
bath a, the evolution equation for the comoving number density Y (= udm/S h) with respect 
to the inverse temperature xi, (= M/Th) of the other thermal bath reads 


dY 

dxb 


{CV) Sb 
XbH 


Y - Y ). 


( 3 ) 


where H is the Hubble parameter, Sb is the entropy density in thermal bath b and l^q 
(= neq/sb) is the equilibrium number density. This equation is not meaningful when Y = 
l^q. This becomes the usual fluid equation in thermal equilibrium. We parameterize the 
annihilation cross section as (av) = cxox^”, in which n = 0 for s-wave annihilation and n = 1 
for p-wave annihilation. The above equation can be reduced to 


dY [tt 


gL 


) - Yl), 


( 4 ) 


where g^, and are the effective relativistic degrees of freedom for the energy density and 
the entropy respectively, and mpj^ is the Planck mass. 

Unfortunately, Eq. (4) has no analytic solution. Fig. 1 shows the result of numerical 
solutions for s-wave annihilation into neutrinos (left panel) and p-wave annihilation into 
photons (right panel). The DM residual annihilation into photons can distort the CMB 


spectrum 


23 


24l |. This effect excludes DMs with mass less than 10 GeV for s-wave an¬ 


nihilation into photons. The effect is negligible for p-wave annihilation, which is velocity 
dependent, so this bound can be evaded. For DM annihilation into neutrinos, we assumed 
that the same numbers of neutrinos and antineutrinos of each type were produced. The 
number of effective relativistic degrees of freedom gl^ is not related to DM or its production. 
We could take the value on the SM base. The p* is taken as a constant g^ on average. The 
curves were made with the proper values of cxop* ' , in agreement with the current DM relic 
density (Iq) for a Dirac fermion, Majorana fermion, complex scalar and real scalar with a 
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FIG. 1; Comoving number density y as a function of inverse temperature x^{= M/T^) and 
Xu{= M/Ty) for a, g = 1 real scalar (short dash), g = 2 Majorana (solid), g = 2 complex scalar 
(dotted) and <7 = 4 Dirac dark matter (long dash) with a DM mass of 10 MeV. The left panel 
is for s-wave annihilation into neutrinos with aog* = (2.8 — 3.0) x 10“^® cm^/s, and the right 
panel is for p-wave annihilation into photons with = 8.0 x 10“^® cm^/s. The horizontal 

dotted line (Tq) represents the current DM relic density, and indicates the equilibrium number 
density. 

DM mass of 10 MeV. As we expected, the DM number track the equilibrium abundance at 
very high temperatures, Xi, < 1. The solution to the Boltzmann equation starts to deviate 
signihcantly from the equilibrium abundance at around x;, ~ 10 — 11, and the actual DM 
abundance Y is different from Veq and Yq for a considerable time. Notice that the equilib¬ 
rium number densities at high tenperatures (early times) are different for different particle 
species. Because the current DM relic density is not relalted to the particle species, a large 
number of relativisitic particles will be produced for the species that has a large equilibrium 
number density at early times. This can be a means to distinguish the nature of the particle 
if other information about the DM, such as the DM mass, is given (see Figs. 2 and 3). 
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B. Thermal Equilibrium Approximation 


DM particles are assumed to keep in thermal contact with one of plasmas after neutrino 
decoupling and to decouple suddenly at some point. The DM and its products can be 
expressed by Fermi-Dirac or Bose-Einstein statistics in this case. There can be different types 
of particles in thermal bath a, so we will use the entropy density Sa = ^g*s'^a = iPa+Pa)/Ta 
to dehne the number of effective relativistic degrees of freedom^ gtsi'^a) in which pa is the 
energy density and pa is the pressure. The energy density pa and the pressure pa are expressed 
by 


Pa = ^Pi = ^ 5^/‘*«5^exp(£)i/r.)±l’ 


(5) 

( 6 ) 


where Pi is the internal degrees of freedom for the corresponding particle i, Ei = \/+ mf 
is the energy with mass and the +(—) sign is for fermions (bosons). We set the chemical 
potentials to zero. The number of effective relativistic degrees of freedom is given by 


T.S (To) 


45 {pa+Pa) 
27r2 Ti 


(7) 


Because the entropy in each thermal bath is conserved after nentrino deconpling, the tem- 
peratnre (Tf,) varies as j, where R is the scale factor. We can find 

the temperatnre ratio at the DM deconpling time if we know the temperatnre ratio at a 
certain time (the time of nentrino deconpling). The temperatnre ratio at the time of DM 
decoupling results in 

Tod _ ( r..iT,,y \ ( gl(T,,o) \ 

no D?.(raD)t VsLWjt ’ 

where is the neutrino decoupling temperature at which the photon and the nentrino 
temperatnres are the same, Tao and are temperatnres^ at DM deconpling. This formula 
can be approximated to the temperatnre ratio at late times; he., at the times when the 
temperatnres T^, Ti, are less than the deconpling temperatnres. Because the DM decoupling 


^ The mark “ ~ ” is placed on top of the symbol of the number of the effective relativistic degrees of freedom 
to indicate DM inclusion. If there is no “ ~ ” mark, DM is excluded. 

^ Notice that one of the DM decoupling temperatures is determined when the equilibrium DM number is 
the same as the present-day DM relic density, Yeq{Ti,D) = T). 
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TABLE I: la and 2a lower limits on the dark-matter mass and upper limits on the existence of 
any other dark radiation for dark matter in thermal equilibrium with neutrinos or electromagnetic 
plasmas. The mark ’ indicates that the limit is irrelevant. The symbol “S” stands for scalar and 
“F” for fermion. 



Neutrino-coupled DM (MeV) 

EM- 

coupled DM (MeV) 


AVeff 


9 

1 (S) 

2 (S,F) 

4(F) 

1(S) 

2 (S,F) 

4(F) 

1 (S) 

2(S) 

2(F) 

4(F) 

la 

6.4 

9.3 

11.9 

10.1 

12.6 

14.9 

0.94 

1.37 

1.27 

1.80 

2a 

3.7 

7.1 

10.0 

- 

- 

- 

1.17 

1.60 

1.50 

2.03 


occurs at ~ 18 as we saw in the subsection A, there must be almost no DM contribution 
to the relativistic degrees of freedom {Tao)- We remove the mark “ ~ . The temperature 
ratio after DM decoupling is then given by 

W.AT.A) ■ 

We now determine Wfr in each case. If the DM particle interacts with a neutrino [a = v 
and 6 = 7 ), the electromagnetic plasma is not relevant to the DM interaction. We can 
identify {g2g{T^) /with {Ty/T-y)^-^. We get the effective number of neutrino 
species from Eqs. (2) and (9) 


SM 


— N' 


9*s {Tyd) 
g'is {Ty) 


4/3 


( 10 ) 


If the DM particle interacts with an eletromagnetic plasma (0 = 7 and b = u), there is only 
one species in the neutrino thermal bath, u. The effective relativistic degrees of freedom 
9*s (Ty) will be the same at any time. Because {Ty/T^)^^ = {gig (T^) /gig we get 

the effective number of neutrino species as 

\9*S D yd) J 


The curves of Fig. 2 display numerical results for the Aefr~ Af relation for a Dirac fermion, 
Majorana fermion, complex scalar and real scalar. The upper (lower) set of curves are for the 
case when DM particles are interacting with neutrinos (electrons and photons). We have 


implicitly assumed that neutrinos decouple at Tyd ~ 2.3 MeV |25l-l27l|. A^"efr increases for 


lighter DM if they are in equilibrium with neutrinos. Conversely, Ws decreases for lighter 
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6.5 


6 

5.5 

5 


g=4, fermion 


g=2, boson 



Dark Matter Mass (MeV) 


10 


20 


FIG. 2: Effective number of neutrino degrees of freedom, as a function of a thermal dark- 
matter mass M. Curves correspond to a (7 = 1 self-conjugate scalar (short dash), g = 2 Majorana 
(solid), g = 2 complex scalar (dotted) and <7 = 4 Dirac dark matter (long dash). The upper (lower) 
curves are for the case when the dark-matter particles are in thermal equilibrium with neutrinos 
(electrons and photons). The dark horizontal band is the Planck CMB Icr allowed range, and the 
light dark band is the 2a upper allowed range. 


DM in equilibrium with an electromagnetic plasma. We put a bound on the DM mass by 
requiring that Wff be compatible with the measured value from Planck j^, and the bounds 
of the DM masses are listed in Table I for each species. If there is a signihcant, but small, 
density of additional radiation, the additional radiation density can be explained by neutrino 
heating from DM annihilation. We should notice that there is still enough room for the 
existence of additional radiation particles (AiVeff), such as a sterile neutrino or a Goldstone 
boson, with decoupling temperatures less than 100 MeV, when DM electromagnetically 
couples to SM particles (EM-coupled DM). 
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C. Out-of-equilibrium Production 


As we can see in Fig. 1, there is a smooth transition between two regions, before and after 
DM freeze-out, and DM particles do not track signihcantly the equilibrium from xi, ~ 10—11. 
The second law of thermodynamics is applied to the entropy (temperature) evolution of the 
produced relativistic particles: 


dSa = 


dQ 

Tl' 


( 12 ) 


where dQ = d{R^p^^) is the heat added per comoving volume due to DM annihilation. 
Because the number of DM particles is reduced by their annihilation at temperature smaller 
than DM mass, the energy density of DM can be described in its nonrelativistic approxima¬ 
tion, ~ — MsbY. The change in entropy^ is given by 


dSa = —SbXbdY —> AS a = —Sb / XbdY 


(13) 


where i is an initial point. We consider the initial point at the time of neutrino decoupling 
because it is the last point at which neutrinos and photons are in thermal contact. DM 
particles with masses less than 20 MeV must be in thermal equilibrium in thermal bath 
a at the initial point. Our observational point is the time of recombination, much later 
after freeze-out. After freeze-out (chemical decoupling), DM continues to scatter off rela¬ 
tivistic SM particles untill DM kinetic decoupling, thus thermalizing the produced particles. 
Produced electrons (positrons) or photons will be thermalized quickly due to the electromag¬ 
netic interaction. Thermalization of neutrinos must be slow. Because neutrinos continue to 
scatter off DM particles after freeze-out, the produced neutrinos^ can be in the equilibrium 
at recombination. The change in entropy is expressed by 


27r2 

AS. = S. - S., = — [gtXR‘ - . (14) 


^ The plasma in the thermal bath b is always in thermal equilibrium because it is not relevant to the DM 
interaction, so the entropy Sb is constant. 

^ A certain number of neutrinos can remain in non-equilibrium if their scattering strength is not enough 
large. We need to consider the detailed Boltzmann equation with the scattering cross section for this 
process. Because our work is not concerned with any specific model, we assume that the produced 
neutrinos are in the equilibrium at recombination. The details with scattering are left for a future study. 
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FIG. 3: Same as Fig. 2, but contour lines are for the case when radiation particles are produced 
from dark-matter annihilation in the non-equilibrium method (freeze-out mechanism). 


The temperature ratio is determined by a combination of Eqs. (13) and (14): 


s) dt)(f) m -ft/--' 


(15) 


where Tai is, according to Ref. jb], very similar to the neutrino decoupling temperature 
described in the SM of the absence of DM. Using the entropy conservation {g'^gR^ ~ T~^) in 
the SM, we can approximate the first term of Eq. (15). The temperature ratio is given by 


rji \ 3 \ 3 

a \ I o, 




'b/ \'^b J SM 9 *s L 

where we have introduced the integration method by parts, as a convenience, for numerical 
computations. The hrst term on the right-hand side of Eq. (16) is just the original tempera¬ 
ture ratio in radiation, and the second term represents a contribution from DM annihilation. 
We can express the temperature ratio in each case, neutrino-coupled DM (a = z/ and 5 = 7 ) 
and EM-coupled DM (a = 7 and b = v). 

Fig 3 shows the numerical results of the iVeff — M relation for neutrino-coupled DM, 
the upper set of curves, and EM-coupled DM, the lower set curves. The basic arguments 
are the same as those in the equilibrium approximation of subsection B. The bounds on 


9ts 


XbY - (xbY). - Ydxb 


(16) 
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TABLE II: Same as Table. I, but the values are for the case when radiation particles are produced 
from dark-matter annihilation in the non-equilibrium method (freeze-out mechanism). 



Neutrino-coupled DM (MeV) 

EM-coupled DM (MeV) 


AVeff 


9 

1(S) 

2 (S,F) 

4(F) 

1(S) 2(S,F) 

4(F) 

1(S) 

2(S) 

2(F) 

4(F) 

la 

9.6 

12.3 

14.8 

9.1 11.8 

14.3 

0.70 

0.99 

0.94 

1.36 

2a 

7.3 

10.3 

12.9 

- - 

- 

0.93 

1.22 

1.16 

1.56 


the DM masses are also listed in Table II, as well as possibilities for the existence of addi¬ 
tional radiation particles (AiVeff). In this DM annihilation process, the DM mass bounds 
are more stringent, and the DM effect on the existence of additional radiation particles 
is stronger. We interpret this in the following way: DM particles annihilate more slowly 
into SM particles to make a smooth transition. In the radiation-dominant era, H = 
{1/R) dR/dt ~ a/(8/3) TiGpn, with the gravitational constant G = l/rrip-^ and the radiation 
energy density p-p . The slower annihilation results in a smaller expansion of the universe and 
eventually a smaller size of the universe later. The same number of relativistic particles must 
be produced from DM annihilation in the equilibrium and the non-equilibrium processes. 
The predicted energy densities at a later time are larger than they are in the equilibrium 
process, so DM annihilation effects are larger in the non-equilibrium process. 


III. CONCLUSIONS 

Light (M < 20 MeV) dark-matter particles freeze out after neutrino decoupling. If 
the dark-matter particle interacts with a neutrino or an electromagnetic plasma, the late¬ 
time entropy production from dark-matter annihilation can change the neutrino-to-photon 
temperature ratio, and equally the effective number of neutrinos Nes- We studied the effects 
of dark-matter annihilation on the N^s by using the thermal equilibrium approximation 
and non-equilibrium method (freeze-out mechanism), and both results were compared with 
Planck observations. If a signihcant, but small, density of additional radiation exists, this 
can be explained by neutrino heating from dark-matter annihilation. The effective number of 
neutrino species Nes is reduced for photon heating. In that case, the existence of additional 
dark radiation particles can help improve the agreement with the current observations. The 
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dark-matter particles annihilate more slowly into SM particles for dark matter annihilation 
in non-equilibrium. The slower annihilation results in a smaller expansion rate (eventually 
a smaller universe later). Although the same number of relativistic particles are produced 
from dark-matter annihilation in the equilibrium approximation and the non-equilibrium 
method, the predicted energy densities at a later time are different. We demonstrated that 
the lower bounds on the dark-matter mass and the possibilities of the existence of additional 
radiation particles are more strongly constrained for dark-matter annihilation process in 
non-equilibrium. 
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